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We derive Tsallis entropy, S q , from universal thermostat independence and obtain the functional 
form of the corresponding generalized entropy-probability relation. Our result for finite thermostats 
interprets thermodynamically the subsystem temperature, Ti, and the index q in terms of the 
temperature, T, entropy, S, and heat capacity, C of the reservoir as Ti = Texp(— S/C) and q = 
1 — 1/C. In the infinite C limit, irrespective to the value of S, the Boltzmann - Gibbs approach 
is fully recovered. We apply this framework for the experimental determination of the original 
temperature of a finite thermostat, T, from the analysis of hadron spectra produced in high energy 
collisions, by analyzing frequently considered simple models of the quark-gluon plasma. 
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A nonlinear entropy formula has been suggested by 
Renyi long ago and been applied to several areas in 
physics [IHZ]. Another formula, the Tsallis entropy, has 
more recently been promoted as the keystone for a gener- 
alized thermodynamics, treating correlated physical sys- 
tems [SHIT]. A respectable amount of papers applying 
this idea to one or the other area in physics appeared 
|12H17j . Since from this entropy the canonical energy dis- 
tribution is power-law tailed in place of the Boltzmann - 
Gibbs exponential, numerous high-energy distributions 
have been fitted using the Tsallis formula |T8H2"T] . Its 
independence from the thermostat and the thermody- 
namical foundation behind the use of such a formula are 
interesting questions. 

In our earlier works we investigated some general 
mathematical properties of alternative entropy formulas 
via their pairwise composition rules, and established that 
a scaled repetition of an arbitrary composition rule leads 
to an associative asymptotic composition rule of large 
subsystems [21]. All such rules are uniquely defined by 
a strict monotonic function, their formal logarithm. Re- 
cently we have also observed that - in connection to the 
zeroth law of the thermodynamics - the factorizability 
condition on the common entropy maximum [29 allows 
only for such rules [30] . We seek in this paper for the 
thermodynamical meaning of the q parameter generaliz- 
ing the classical entropy formula, valid for q = 1. Some 
q 7^ 1 parameter was calculated theoretically [3TH55] . 

We found the thermodynamical interpretation of the 
entropy formula and its parameters on the analysis of 
the two-body thermodynamics of a single observed sub- 
system and a reservoir. For finite systems the micro- 
canonical approach is the key to the physical interpre- 
tation. In the classical treatment subleading terms in 
a finite-energy expansion of the microcanonical entropy 
maximum are usually ignored, the reservoir is treated as 
constant in the canonical limit. A notable exception is 
the analysis of statistical fluctuations and their scaling in 
the thermodynamical limit [34 440j . We consider the cor- 
relation between subsystem and reservoir induced by the 



conservation of total energy while maximizing a mono- 
tonic function of the Boltzmann - Gibbs entropy, L(S). 
We seek for that very function, L, which counteracts fi- 
nite size effects beyond the usual linear term, —j3Ei, in 
the Taylor-expansion of the L(S) — max principle. 

We discuss now the thermal equilibrium of two sys- 
tems, one with energy E\ (subsystem) and the other with 
energy E — E\ (reservoir) , while their respective entropy 
contributions are combined by the general rule satisfying 

L(S 12 )=L(S 1 )+L(S 2 ). (1) 

Here we do not assume that the deviation from the sim- 
ple additive rule would be small. For the sake of simplic- 
ity we consider here homogeneous rules, relevant for the 
cases when subsystem and reservoir are composed from 
the same matter (for details see Ref. 30J ) . The micro- 
canonical condition for a maximal entropy state then de- 
fines the thermodynamical inverse temperature, requir- 
ing 

L(S(E 1 )) + L(S(E-E 1 ))= max. (2) 

Varying the subsystem energy, E±, while keeping the 
total energy E fixed, we describe the thermal contact 
between subsystem and reservoir. This means that the 
derivative with respect to E\ of the above expression ^ 
vanishes. Owing to the two -Ei-dependent contributions, 
it is equivalent to the statement that 

ft = L'(S(E 1 ))-S'(E 1 ) = L'(S(E-E l ))-S'(E-E 1 ). (3) 

This equality, when taken in the E 3> Ei limit, usually 
defines the canonical approach. Now we would like to 
take into account effects to higher order in E\/E, and 
require that their leading term vanishes on the right hand 
side. The reservoir's entropy on the right hand side is 
Taylor-expanded: S(E - E x ) = S{E) - S'(E) ■E l + . . .. 
Collecting the coefficients of E\ we arrive at 

ft = L'(S(E)) ■ S'(E) 

- [S'(E) 2 L"(S(E)) + S"(E)L'(S(E))] E 1 + ... (4) 
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Here the first term on the right hand side is the familiar 
canonical (E'l-independent) Lagrange multiplier, 



/3 = L'(S(E)) ■ S'(E) 



L'(S) ■ \ 



(5) 



constituting the /3i = j3 relation. Our key addition to the 
usual treatment is to require that the coefficient of the 
linear term in eq.Q vanishes: This is a constraint for the 
L(S) function in general. Obviously, without considering 
L(S), the whole coefficient consisted only of S"(E) as in 
the traditional approach, and nothing further could be 
done. We obtain the following condition: 



L"{S) _ S"(E) 
L'(S) ~ S'{E) 2 



(6) 



Since the left hand side of eq.Q is a function of S, 
while the right hand side is a function of E, the left hand 
side must be treated as an jS-independent constant by 
solving eq.([6]) for L(S). This Universal Thermostat In- 
dependence (UTI) reads as 



L"(S) 
L'(S) 



(7) 



The particular solution with L'(0) = 1 and L(0) = is 
given by 



L(S) = 



e aS_ l 



(8) 



The derivatives of the S(E) equation of state do have 
physical meaning: S'(E) = l/T and S"{E) = -l/CT 2 
are related to the traditional temperature and heat ca- 
pacity of the reservoir. By using this we obtain 



a = 1/C. 



(9) 



The non-additivity parameter is simply the inverse heat 
capacity of the reservoir. For C — > oo one has a — > and 
L(S) — > S, so the Boltzmann - Gibbs formula is included 
by this limit. A connection between Tsallis entropy and 
constant heat capacity of the reservoir has been observed 
years ago [HJ ■ Our result shows the background for 
this observation. The philosophy behind our approach is 
first to decide on the entropy formula by choosing L(S) 
generally and then to solve the maximization problem in 
terms of subsystem energies and corresponding probabil- 
ities. 

The knowledge gained from this analysis now will be 
generalized. Analog to a Gibbs ensemble, we extend a 
sum of two to a weighted sum of many. In this way the 
result of the two-body analysis in the form L(Si) gen- 
eralizes the classical entropy formula, S = — J2i Pi m -Pj> 
to 



L(S) = Y,PM-\nPi). 



(10) 



This L- additive form of the generally non- additive en- 
tropy leads to the Tsallis entropy formula when applying 
eq.([8|. In this way one obtains 

L (S(Ei)) -f3E 1 = - (e aSi - E ^ - l) - 0Ei = max. (11) 



Here the coefficient of the second order correction, 
O (El/E 2 ), vanishes for a = l/C{E). By this we are 
led to the following entropy expression 

L(S(E 1 )) = L(-lnPx) = - a (Pr a - 1) ■ (12) 

Our result applied to a Gibbs-ensemble with the relative 
occurrence frequency Pi of states with energy Ei, hence 
reads as 

^ PM- In Pi) - p P i E i ~ a J2 Pi = max ' ( 13 ) 



Substituting eq.(12) we finally arrive at 



I Yl ( p *~ a - p )-pY. ~ a J2 p * = max - ( 14 ) 

i i i 

With the widespread notation q = 1 — a one obtains the 
Tsallis entropy formula 



tsallis :=L{S) = — i Y,{P*-P?)- 



(15) 



It is suggestive to consider its inverse function, L 1 ac- 



cording to eq.(10). This delivers the Renyi entropy: 

l-o ^ 1 



S 



Renyi 



s 



(16) 



Now the parameters j3 and a, defined in eqs. @ and j9j, 
are set by the physics of the finite-energy reservoir. The 
sign of the heat capacity, C, determines whether the q 
is smaller or larger than one. It may possibly carry an 
interesting message for the description of gravitating sys- 
tems, with negative heat capacity. 

We proceed by noting that maximizing S'xsaiiis with 
respect to the Pi weights of system instances with energy 
Ei one obtains the canonical cut power-law distribution 
of energies: 



P,= I Z 1 -" + (1 - q)^Et 



(17) 



Using eqs. ([5]), and ([9| we rewrite this in the equiva- 
lent form, 



Pl ~z( 1+ C-l T 



(18) 



expressing the energy distribution in terms of the tem- 
perature, T, entropy, S and heat capacity, C of the ideal 
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reservoir. The partition sum Z, obtained from normaliza- 
tion, is related to the Tsallis-entropy, L(Si), and energy, 
Ei, of the subsystem via its deformed logarithm: 

In, Z := C - l) = LiS,) - —L^l (19) 

In the infinite heat capacity limit, irrespective to the 



of the energy with respect to temperature, 



value of 5, formula ( 18 1 recovers the exponential dis- 



tribution. The inverse logarithmic slope of the energy 
distribution, derived from it, is linear: 



/ / ''"d/v 1 "''*' =Tu + E,;C. (20) 



with T = Te- s l c Z 1 l c {\ - 1/C). One concludes that 
the generalized entropy formula leads to a cut power-law 
energy distribution, based on a finite heat capacity reser- 
voir. As such, it is a better approximation to the micro- 
canonical distribution, than the canonical exponential. 

We demonstrate the usefulness of the above general re- 
sults on the example of the thermal model to heavy-ion 
collisions. Experimental data from RHIC AuAu colli- 
sion at 200 GeV deliver different T slopc -s extrapolated 
to px — for different hadrons [THl EH|- Considering 
that the energy at zero momentum is the rest mass in 
c = 1 units, a linear trend shows in the T s i opc (mi) val- 
ues, as seen in Fig[l] The open circles correspond to 
mesons, the filled ones to baryons in this figure. The 
steepness for mesons and baryons seem to be in the pro- 
portion 2 : 3 suggesting a quark coalescence hadroniza- 
tion picture, compatible with the factorization assump- 
tion P h adron(£) = P* ' (E / ' K) with K = 2 and K = 3 
for mesons and baryons, respectively. This scaling is ac- 
ceptable and leads to T^° n {E) = T^ c k (E/K), with 
the common To as 48 MeV intersect in the formula ( 20 1 . 



The valence quark matter heat capacity at RHIC AuAu 
collision tends to be C ~ 4.5. Similar trends can be ex- 
tracted from the analysis of fits to the ALICE data in 
900 GeV pp collisions done in [?3]: the values for the 
Tsallis-slope parameters, To, are much lower than the 
canonical QCD phase transition temperature. Here we 
re-plotted the tabulated values given in [44] using the 
coalescence quark assumption, denoting mesons by open 
square boxes while baryons by filled boxes. We note, 
however, that these fits were performed in the very low 
Pt range ( pr < 2.5 GeV/c) only, therefore the uncer- 
tainty of the fitted parameters is large. 

In order to interpret this surprisingly low value for 
To, we have to consider physical models of a finite ther- 
mostat, and calculate T\ = = Te~ s l° ', since 
lim T) = Ti for small subsystems in large reservoirs. 

C— foo 

First we study the Stefan - Boltzmann formula supple- 
mented with a bag constant, E/V — aT 4 + B in a volume 
V. Since the pressure is given by p — \oT 4 — B, the en- 
tropy is 5 = iaVT 3 . The heat capacity is the derivative 



C= ^ =W T 3 
dT 



(<xT 4 + B 



dV 
d~T' 



(21) 



At constant volume, V, this gives Cy = iaVT 3 = 35 
and T±v = Te -1 / 3 . At constant pressure the temper- 
ature cannot change in this model, so C v — oo and 
Tip = T. Furthermore, considering an adiabatically ex- 
panding reservoir, a more realistic scenario in high-energy 
experiments, one deals with the heat capacity at con- 
stant entropy, C s = 35(1 - T 4 /T 4 )/4, with T being the 
temperature where the pressure vanishes. In this case 
Cs < 35/4 and T 1S < Te~ 4 / 3 is the theoretical predic- 
tion. 

Figure [l] presents the inverse logarithmic slope, T s i opc , 
as a function of hadron masses and Ti-lines for differ- 
ent physical models of the thermostat. Besides the three 
above described bag-model approaches we also indicate 
the classical Schwarzschild black hole, having C = —25 
and Ti = Te 1 / 2 . One inspects that this possibility is far 
from all experimental observations. 

We note that theoretically a really constant heat ca- 
pacity, Co, stems from the equation of state S(E) — 
Co ln(l + c^Tb )■ The latter is a good ansatz for an effec- 
tive equation of state of classical non-abelian gauge field 
systems on the lattice [45] and represents the high-T limit 
of Planck's 5"(T)-formula for thermal radiation. 

Considering heat capacity in the above scenarios and a 
standard numerical value of T « 167 MeV for the reser- 
voir temperature, conjectured for the QGP at hadroniza- 
tion phenomenology and determined by lattice QCD cal- 
culations, one obtains T 1P = T = 167 MeV, T 1V = 
Te" 1 / 3 a 120 MeV and T 1S < Te~ 4 / 3 w 45 MeV char- 
acterizing the Tsallis-distribution of valence quarks. 

The conjecture that in heavy ion collisions a statis- 
tical power-law energy distribution due to finite phase 
space availability corrections to the traditional canonical 
distribution may appear is further supported by the ob- 
servation that the measure of non-additivity, a = 1/C, 
expressed by the inverse power in the fitted power-law 
tail, is reduced for increasing participant number |46) . 
The fitted power C is also tendentiously smaller in e + e~ 
or pp than in heavy ion collisions [21] . Finally we real- 
ize that only the adiabatic scenario for the quark matter 
thermostat leads to Tq values near to the ones extracted 
from experimental analysis by the coalescence assump- 
tion, T 1S a T a 45 - 55 MeV. 

In conclusion the Tsallis entropy formula is derived as 
the consequence of the following requirement: we seek for 
that non-additive entropy composition rule which cancels 
linearly energy-dependent corrections due to the finite 
E — Ei energy in the reservoir to a subsystem's ther- 
modynamical inverse temperature. This determines the 
composition rule and the entropy formula uniquely turns 
out to be the Tsallis entropy. This derivation explains the 
particular functional form of the Tsallis and Renyi for- 
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FIG. 1. Extrapolated inverse slopes from RHIC AuAu data 
at 200 GeV (circles), LHC pp data at 900 GeV (boxes) and 
for different theoretical models of a QGP-thermostat with 
temperature T = 167 MeV (horizontal lines). Open sym- 
bols correspond to mesons, filled symbols to baryons. The 
skew lines indicate the respective valence quark assumptions: 
TZ;T(pt = 0) = To + mi /2C and T^TiPT = 0) = 
To + rrii/3C, for the LHC (upper two lines) and for RHIC 
(lower two lines). 



mulas as generalized entropy expressions satisfying the 
UTI principle. With regard to the physical interpreta- 
tion we have obtained q = 1 — 1/C, with C being the 
heat capacity of the total system with the conserved en- 
ergy E. The canonical temperature of the subsystem 
becomes 7i = e~ s /° T with T(E) and S(E) being the 
traditional temperature and entropy of the finite reser- 
voir, respectively. A preliminary analysis of experimental 
data on particle production seems to be sensitive to dif- 
ferent physical assumptions about a QGP thermostat. 
Here the iscntropic scenario performs best. 
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